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The quasiparticle interferences (QPIs) of the featureless Mott insulators are investigated by a 
T-matrix formalism implemented with the dynamical mean-held theory (T-DMFT). In the Mott 
insulating state, due to the singularity at zero frequency in the real part of the electron self energy 
(ReE(a;) ~ p/oj) predicted by DMFT, where r) can be considered as the ‘order parameter’ for the 
Mott insulating state, QPIs are completely washed out at small bias voltages. However, a further 
analysis shows that ReE(a;) serves as an energy-dependent chemical potential shift. As a result, 
the effective bias voltage seen by the system is eV' = eV — ReE(ey), which leads to a critical bias 
voltage eI4 ~ ,Jr\ satisfying eV = 0 if and only if t] is non-zero. Consequently, the same QPI 
patterns produced by the non-interacting Fermi surfaces appears at this critical bias voltage eI4 
in the Mott insulating state. We propose that this re-entry of non-interacting QPI patterns at eVc 
could serve as an experimental signature of the Mott insulating state, and the ’order parameter’ can 
be experimentally measured as ~ (eI4)^. 

PACS numbers: 71.27.-|-a,72.10.Fk,71.10.Fd 


Introduction - Mottness, the physics of understanding 
how the insulating state arises from a partially-filled con¬ 
duction band due to strong local interactions, has been 
one of the most challenging subjects in condensed mat¬ 
ter physics. It could potentially hold the key to under¬ 
standing the mechanism of high-temperature supercon¬ 
ductivity observed in various materials with a narrow 
bandwidthi^— . The essential physics of Mottness can 
emerge from the single orbital Hubbard model, 

^Hubbard — Y Hjj 

Ht = -t 4aCj<T + h.C., 

<i,j>,<T 

Hu = (1) 

i 

where t is the hopping parameter between nearest neigh¬ 
bor sites < i,f >, and U is the on-site Coulomb interac¬ 
tion. Although we analyze a single orbital model here, 
our conclusions are also applicable to multi-orbitals sys¬ 
tems. While the Hubbard model has been successful in 
conceptually demonstrating both the metallic and the 
insulating states in two extreme limits oiUjt « 1 and 
U/t >> 1, the nature of the transition from the metallic 
to the insulating states as a function of U/t is still ob¬ 
ject of active investigation. One particular difficulty is 
that the Mott transition does not seem to involve break¬ 
ing of symmetry, and consequently such a transition can 
not be studied within the framework of Landau-Ginzburg 
theory. The change in the temperature dependence of 
dc resitivityi is often used as an indicator of the Mott 
transition, but it depends on many other details of the 
materials, e.g., impurity and disorder. As a result, the 
dc resitivity can not be a conclusive experimental signa¬ 
ture for the Mott transition. Other possible candidates 
for an ‘order parameter’ of the Mott transition, such as, 


the density of states at Fermi energy^ or the number 
of double occupancy^, have been discussed in the liter¬ 
ature. However, these physical quantities are also not 
ideal candidates for an order parameter behavior since 
for example, the density of states at Fermi energy is zero 
in the insulating state and therefore does not contain in¬ 
formation about the strength of Mott insulator. In other 
words, an order parameter for the Mott transition that is 
theoretically well-defined and experimentally observable 
is lacking. 

The dynamical mean-field theory (DMFT) has been 
shown to be a sophisticated theoretical approach to un¬ 
derstand the Mott insulating state^Ai. DMFT maps the 
problem of solving the Hubbard model to another one 
of solving the Anderson impurity model which is much 
better understood. Such a mapping is shown to be ‘ex¬ 
act’ in the limit of spatial dimension d going to infin¬ 
ity. The price paid for this mapping, however, is that 
the momentum-dependence of the electron self energy is 
completely neglected, but the resulting electron self en¬ 
ergy E(a;) includes all the local quantum fluctuations. 
Since the Mott insulating state concerns mostly the local 
physics, E(a;) from DMFT captures a lot of information 
that can not be accessed by other approaches easily, de¬ 
spite its lacking of the momentum dependence. 

In DMFT for the case of half-fillingi, the real part of 
the electron self energy has the following singular behav¬ 
ior at low frequencyi^ 

ReE(..) = ^+0(cc), (2) 

Here rj is given by, 

/ OO _ 

^^ImGb(fc,e-biO+)] \ (3) 
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where Gb{k,ui) is the Green function obtained from 
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DMFT with the following form, 

Gb{k,uj)=[uj-E{k)-i:iuj)]~\ (4) 

E{k) = e{k) — where e(fc) = —2t{coskx + cosky), and 
/i is the chemical potential. The constant term —U/2 has 
been absorbed into the chemical potential. 

This l/w singularity in ReE(w) is inherited from the 
self energy in the atomic limit of the Hubbard model 
(t = 0), and DMFT captures it well even with non-zero 
t. It can be shown^ that r; = 0 in the metallic phase while 
ry ^ 0 in the Mott insulating phase. Moreover, a larger 
T] corresponds to a stronger Mott insulating state. Due 
to the above reasons, rj has been proposed to have prop¬ 
erties similar to an order parameter of the Mott insulat¬ 
ing stateiii^. However, r] remains a purely theoretically- 
defined ‘order parameter’ since no experimental probe 
has been applied to measure it. 

In this paper, we propose that the quasiparticle in¬ 
terferences (QPI’s) from impurities on a material sur¬ 
faces which are measured by spectroscopic imaging scan¬ 
ning tunelling microscopy (SI-STM) could be a feasi¬ 
ble tool to directly measure rj. QPIs have been widely 
used to extract the spectroscopic information of the 
electronic structure in materials including, but not lim¬ 
ited to, cuprateai^^— , iron-based superconductors'^^— , 
topological insulators^Sr— , Sr 3 Ru 2072 i“^, and heavy 
fermion systems^ii^. On the theoretical side, the T- 
matrix formalis m^^i^^i^^i^° has been shown to be a well- 
established method to compute the QPI images. We de¬ 
velop a T-matrix formalism implemented with DMFT 
(T-DMFT) to analyze the QPI patterns in the Mott in¬ 
sulating state. We find that at small bias voltage, no 
meaningful QPI images are found due to the singular 
behavior of the electron self energy given in Eq. [2l as ex¬ 
pected. However, by analyzing the T-DMFT formalism, 
we further find that the QPI patterns resembling the one 
produced by the non-interacting Fermi surface appear at 
a non-zero critical bias voltage Vc, and this critical bias 
voltage is directly related to ij, the ‘order parameter’ of 
the Mott insulating state given in Eq. [51 We propose 
that this novel re-entry of non-interacting QPI at high 
bias voltage could serve as an experimental signature of 
the Mott transition, and the order parameter can be ob¬ 
tained by 77 ^ (eVc)^ accordingly. 

Formalism - We start from the following model Hamil¬ 
tonian of 

ki E}jxibhard T 

Himp = X/X/ 

k,k' ^ 

where Hnubbard is the Hubbard model described in Eq. [1] 
and Vg p is the impurity scattering matrix element which 
is assumed to be spin-independent. To specifically com¬ 
pute the QPI image, we employ a T-matrix approach22i2£ 
implemented with DMFT, and the procedure is given be¬ 
low. 


Following Ref. Q , we first solve H Hubbard with DMFT 
using Hirsch-Fye quantum Monte Carlo algorithm as the 
impurity solver. The imaginary time interval [0, /3] is di¬ 
vided into L = 128 slices, and the maximum Matsubara 
frequency used is Umax = 8192. We use the semicircular 
density of states with half-bandwidth D. After the self- 
consistent calculation is converged, the time-ordered elec¬ 
tron self energy in Matsubara frequency S(fa;„) can be 
evaulated directly. In order to obtain the electron Green 
function in real frequency Gb{k,uj), a numerical calcu¬ 
lation of the analytic continuation on is neces¬ 

sary. We adopt the continuous-pole-expansion method^! 
recently developed by Staar, et al. to numerically per¬ 
form the analytic continuation on This method is 

shown to be an efficient approach to find an accurate and 
unambiguous result of the analytic continuation within a 
finite range of frequency, which is necessary for the cur¬ 
rent study. Finally, the electron Green function with 
DMFT Gb(k,uj) given in Eq. Ujcan be obtained after the 
analytic continuation on E(za;„) is performed. 

Once Hnubbard is solved by DMET, we treat impurity 
scatterings necessary for the QPI images as a perturba¬ 
tion around the DMET solution. If the impurity scatter¬ 
ings are weak, T-matrix formalism22i^ allows us to sum 
over all the vertex corrections due to the impurity scat¬ 
terings with Gb{k,uj) as the ‘bare’ Green function. The 
full Green function including impurity scatterings can be 
expressed 

G{k,k',io) = Gb{k,uj)6^p F Gb{k,uj)T^p{u:)Gb{k',uj) 

( 6 ) 

where the T-matrix satisfies 

+ Jd^pVj: .Gbip,u;)T.j:,{u:) , ( 7 ) 

Typically, SI-STM is performed on a nicely-cleaved sur¬ 
face on which the impurities are usually apart from each 
other. In this case, the STM conductances (^) at po¬ 
sitions near a single impurity are measured, and conse¬ 
quently it is sufficient to include only one impurity in 
the calculatio n^^’^° . Eor the case of the single impurity, 
VjT p = Vo, and can be further reduced to T(a;). 

As a result, T(w) can be evaluated from Gb{k,uj) by 

= ^ - J(fk Gb{k,ui) (8) 

It is well-establishe d^^i^° that ^{f,eV) at bias volt¬ 
age V is proportional to the local density of states 
plr,eV) = where Pa{r,uj) = lmGa{r,r,uj). 

The QPI image in the q map , namely p{q, w), is just the 
Eourier component of p{r,eV). After a straightforward 
calculation, the T-DMET formalism leads to 

piq,eV) - PoieV) = 

lm{T{eV) Jd^kGbik,eV)Gbik + q,eV)}, ( 9 ) 
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FIG. 1: The Fermi surface of Hnubbard with t/ = 0 at half¬ 
filling. The red dashed lines represent the Fermi surfaces, and 
the blue arrow indicates the momentum (in this case, (7r,7r)) 
connecting two points on the Fermi surface with large joint 
density of states. 
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where po{eV) is the density of states without the impu¬ 
rity scatterings which will be dropped out from now on 
since it is just an unimportant constant. 

Results - For the non-interacting case, S(w) in Eq. 0] 
is zero, and the peaks in the q map of p{q, eV) are cor¬ 
responding to the momenta connecting two points with 
large joint density of states on the equal energy contour 
at energy eV away from the Fermi energy22i2£. For ex¬ 
ample, the Fermi surface of Hnubbard with [/ = 0 at 
half-filling is depicted by the red dashed lines in Fig. [T] 
The q map of p{q, eV = 0) has a peak at q = (tt, tt) 
corresponding to the blue arrow in Fig. [TJ 
As the interaction is present, 

Gb{k,uj) = [uj-E{k)-j:{oj)]~^ 

= [w'— if(fc) — zImS(a;)] (10) 

where oj' = uj — ReE(eE). Consequently, Eq. IHlbecomes 
piq,eV)=lm{T{eV) (11) 

_ ££ _I 

[eV-E(k)-ilmY.{eV)][eV'-E{k+^-ilmY.(eV)]k' 

It is now clear that eV = eV — ReS(eE) is the effective 
bias voltage shifted due to the present of ReE(eR). The 
QPI images at bias voltage eV in this case is in fact the 
same as the non-interacting one at eV' except the extra 
broadening resulted from ImE(eR). 

The electron self energy E(w) obtained from DMFT 
calculation for various valules of U is plotted in Fig. [2] 
Note that both U and ksT are in the unit of D as used 
in Ref. 0- For small C/, the self energy is small, and 
therefore the QPIs are not expected to deviate from the 


FIG. 2: The real (top) and imaginary (bottom) parts of the 
electron self energy from DMFT for U = 1,2, 3, 4, 5 with tem¬ 
perature ksT = 0.1. Note that the unit of energies is D as 
used in Ref. 0. The critical bias Vc is the point crossed by 
the line of uj and ReE(aj). For 17 = 1, 2, the system is not in 
the Mott insulating state, thus Vc = 0. 

ones obtained in the non-interacting case. The situation, 
however, changes dramatically as U is large enough for 
the occurrence of the Mott insulating state (2 < C/c < 3 
in the current calculation). In the Mott insulating state, 
ImE(a;) has a delta-function-like peak at u = 0, while 
ReS(a;) exhibits a divergence of 1/w at small frequency. 
These behaviors of self energy make Gb{k,uj) very small 
at small frequency, which washes away any meaningful 
features in the QPI images at small bias voltage. 

Intriguing features appear at high bias voltage. The 
singular behavior of ReE(a;) described in Eq. 0] guaran¬ 
tees the existence of a critical bias voltage eI 4 satisfying 

eVf = et4 - ReE(eI4) = 0. (12) 

Adopting the expansion of ReE(a;) given in Eq. [21 we find 
eRc ~ y/p. Although the precise value of eI4 depends on 
the high order corrections, it is important to recognize 
that eVc scales with ij, and eI 4 = 0 if and only if 77 = 
0. In other words, in the Mott insulating state, there is 
always a non-zero critical bias voltage 14 at which the 
QPI pattern is exactly the same as the non-interacting 
one except extra broadening, and the ‘order parameter’ of 
the Mott insulating state can be indicated by rj (el4)^. 

To confirm this, we consider the one-band Hubbard 
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model with the nearest neighbor hopping at half-filling. 
The non-interacting QPI exhibits a strong peak at g = 
(tt, tt) at zero bias and the peak shifts to a smaller mo¬ 
mentum as the bias voltage increases. Fig. [3] plots the 
QPI patterns along the nodal direction in the first Bril- 
louin zone at the critical bias voltage 14 for different U. 
All of them exhibit exactly the same peak at q = (tt, tt), 
resembling the non-interacting QPI at zero bias. For 
U = 1,2, the system is not in the Mott insulating state 
and 14 = 0. For U = 3,4,5, the non-interacting QPI 
pattern appears at el4 = 0.44,0.88,1.37 respectively. 
Longer range hopping parameters can be added in or¬ 
der to capture more details in the Fermi surface, but 
these additions do not change the conclusions except the 
critical bias voltage might be slightly modified. 

Discussions and summary - The current result is pri¬ 
marily based on the feature that the electron self energy 
obtained from DMFT is momentum-independent. This 
feature is correct in the atomic limit of the Hubbard 
modeli, but in general the electron self energy should 
depend on the momentum whenever t/U is non-zero. 
Nevertheless, since we are mainly interested in the Mott 
insulating state in which t/U is a small parameter, the 
momentum-dependence of the electron self energy should 
be very weak in the Mott insulating state, which has been 
discussed in previous studies^^— . As a result, the pre¬ 
dicted re-entry of the non-interacting QPIs at a critical 
bias voltage should be robust even as the momentum de¬ 
pendence of the electron self energy is considered. 

It is worth being mentioned that the critical bias volt¬ 
age is smaller than the Mott gap as we can see in Fig. |31 
The typical value of D for cuprates is around 400 meV, 
which gives the critical bias voltage about 150 meV to 600 
meV for U = 3 — 5. The fundamental limitation on the 
energy resolution of SI-STM is the range of the energy 
within which the density of states of the STM tips re¬ 
mains flat. A typical STM tip has a flat density of states 
within Ep ± 500 meV. As a result, although the critical 
voltage is a little higher than the typical range of the 
bias voltages where the SI-STM is routinely exploited, it 
is still measurable with the current techniques. 

The existence of the critical bias voltage eI4 is not 
limited to the electron self energy obtained from DMFT. 
In fact, any theory having an electron self energy with 
a diverging behavior at low frequency is likely to have a 
non-zero eI4- From the Dyson’s equation, the self energy 
can be expressed as E = Gq^ — G“^, where Go and G are 
the non-interacting and the full Green functions. Since 
Go is non-zero always, E could have a diverging behavior 
if and only if the full Green function has zeros at low fre¬ 
quency. The existence of zeros in full Green function has 
been recently found in a variety of theories for strongly 
correlated electron o^^i^^~ — , including the AdS-CFT holo¬ 
graphic models^i^. As a result, it is expected that the 


same result will be obtained if these theories could be 
incorporated into T-matrix formalism to compute QPIs. 
Further, note that although for computational simplicity 
we have focussed on the 2D Hubbard model which is ap- 

QPI at eVc 
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FIG. 3: The QPI at the critical bias voltage Vc for dif¬ 
ferent U along the nodal direction. p{q, eI4) is normalized 
with respect to the peak value at g = (tt, tt) for each U. 
eI4 = 0, 0, 0.44, 0.88,1.37 for U = 1,2, 3, 4, 5. 


plicable to a large class of materials including cuprates 
and Heavy-Fermions, our conclusions are more general 
and applicable to 3D systems as well. 

In summary, we have employed the T-matrix formal¬ 
ism implemented with the dynamical mean-field theory 
(DMFT) to study the QPI patterns of the featureless 
Mott insulator. While QPIs at small bias voltages are 
completely washed out due to the singular electron self 
energy obtained from DMFT, we find that the QPI pat¬ 
terns resembling the non-interacting ones appear at a 
non-zero critical bias voltage. Since the existence of this 
non-zero critical bias voltage is a direct consequence of 
the singular behavior of ReE(a;) ^ ij/uj in which rj could 
be thought of the ’order parameter’ of the Mott insu¬ 
lating state, this novel re-entry of non-interacting QPI 
patterns at a finite critical bias voltage could serve as an 
experimental signature of the Mott transition, and the 
order parameter can be indicated by 77 ~ accord¬ 

ingly. 
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